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^ ' Abstract 

^^ . 

\^ , The intermittent transition between slow growth and rapid shrink- 

lO ' age in polymeric assemblies is termed dynamic instability, a feature 

^*l \ observed in a variety of biochemically distinct assemblies including 

t^^ ■ microtubules, actin and their bacterial analogs. The existence of this 

^D I labile phase of a polymer has many functional consequences in cy- 

toskeletal dynamics, and its repeated appearance suggests that it is 
relatively easy to evolve. Here, we consider the minimal ingredients 
for the existence of dynamic instability by considering a single poly- 
morphic filament that grows by binding to a substrate, undergoes a 
''^ • conformation change, and may unbind as a consequence of the resid- 

ed I ual strains induced by this change. We identify two parameters that 

control the phase space of possibilities for the filament: a structural me- 
chanical parameter that characterizes the ratio of the bond strengths 
along the filament to those with the substrate (or equivalently the 
ratio of longitudinal to lateral interactions in an assembly), and a ki- 
netic parameter that characterizes the ratio of time scales for growth 
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and conformation change. In the deterministic hmit, these parame- 
ters serve to demarcate a region of uninterrupted growth from that of 
coUapse. However, in the presence of disorder in either the structural 
or the kinetic parameter the growth and cohapse phases can coexist 
where the filament can grow slowly, shrink rapidly, and transition be- 
tween these phases, thus exhibiting dynamic instability. We exhibit 
the window for the existence of dynamic instability in a phase diagram 
that allows us to quantify the evolvability of this labile phase. 

Key words: mechanochemistry ; protein polymorphism ; dynamic 
instability ; microtubule ; actin ; elastic filament 

Introduction 

Polymeric filaments are the building blocks of nearly all biological 
structures at the cellular level. Their structural, chemical and mechan- 
ical properties control several processes in the cell and beyond. A key 
property of these filaments that allows them to be so flexible in their 
structure is their dynamic lability which allows them to grow or shrink, 
and become crosslinked or fall apart with relative ease. This is achieved 
through a variety of structural and chemical means such as capping, 
treadmilling and most spectacularly, dynamic instability. This last 
process which involves intermittent transitions between phases of slow 
rowth and rapid shrinkage of a polymer, first seen in microtubules 

After the discovery of dynamic instability in microtubules llj), var- 
ious models ^]m) have been proposed to interpret the phenomenon as 
a stochastic process with different kinetic constants for the addition 
and removal of subunits from the ends of polar subunits. While this 
chemical kinetic approach leads to results that are able to explain the 
experimental observations qualitatively, over the years it has become 
increasingly clear that dynamic instability in microtubules, the best 
studied system to date, has an important structural component asso- 
ciated with the change in shape of the dimers once the attached GTP 
is hydrolyzed (|6|-l8|). In particular, microtubules which are formed by 
a number of protofilaments, typically thirteen, grow by the addition 
of tubulin dimers which are in their GTP-bound state. Soon after 
polymerization, the GTP-bound tubulin changes conformation from a 
straight state to a bent state upon the hydrolysis of the GTP unit 
([y, |8|). This conformational change is critical to dynamic instabil- 
ity, since curved filaments tend to detach from the microtubule while 
straight filaments are stable. Indeed electron micrographs of micro- 
tubules caught "in flagrante delicto" show that individual protofila- 
ments can be seen curving outwards from the frayed ends (121), and 
more recently individual protofilaments have been found to assemble 
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into rings that curve along a direction orthogonal to that when they 
are part of the tubule (|8J) . Mechanical measurements of the rigidity of 
microtubles Q) show that the Young's modulus of the assembly is two 
orders of magnitude smaller than its shear modulus, consistent with 
the structural evidence of strong interactions between tubulin dimers 
along a protofllament and weak lateral interactions between tubulin 
dimers on different protofilaments. Taken together these observations 
suggest that the stability and dynamic instability of microtubules in- 
volves structural, mechanical and kinetic aspects (|l0l - [l4l ). Large scale 



computationally intensive models (jlSI ) do try and to account for these 
effects, but at the cost of understanding the generality and qualitative 
aspects of the basic phenomenon. 

Since the observations of dynamic instability in microtubules more 
than 20 years ago, the phenomenon has been implicated in the dynam- 
ics of single actin filaments (jlGf ) , observed in the bacterial homolog of 
actin, ParM (1171), and is thought to also occur in bacterial homologs of 
microtubules (|l8l ). In all these cases, the process of dynamic instabil- 
ity is characterized by the cooperative assembly and disassembly of a 
filamentous multi-stranded polymer that involves multiple structural 
states of the subunits and filament ends (l3) that are dependent on the 
hydrolysis of GTP or ATP. Thus, the lability of the polymeric filament 
is directly tied in to the kinetics of assembly as well as the consequent 
structural polymorphism, and conveniently used by the cell to direct 
its internal organization during movement, division and other func- 
tions. The appearance of similar kinetic profiles for the growth and 
shrinkage of polymeric filaments that are very different biochemically 
suggests a functionally-driven convergent evolution that selected cer- 
tain traits in subunits capable of self assembly. This naturally raises 
the question of the minimal system that is capable of robust dynamic 
instability in polymeric assemblies with the following ingredients (i) 
a slender geometry associated with both the individual protofilament 
and the filament assembly, (ii) bond interactions of different strengths 
along and across protofilaments, (iii) kinetics of subunit addition being 
different from that of NTP hydrolysis in the added subunit. Each of 
these requirements is fairly generic and suggests a model of a single 
elastic filament which can grow or shrink by subunit addition, can at- 
tach to or detach from a substrate (of other protofilaments) and has in 
internal degree of freedom associated with an NTP-hydrolysis driven 
conformation change. Our model is similar in spirit to the filament 



model introduced in Ref. (jllI ) 



The relative simplicity of the model allows us to use a combination 
of scaling arguments, analytic solutions, and simple numerical simula- 
tions to characterize both the deterministic and stochastic aspects of 
the growing and shrinking phases of the elastic filament. In particular, 
we find that for the completely deterministic case, simple arguments 
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allow us to characterize the mechanochemical conditions under which 
filaments either grow without bound or collapse and disappear. How- 
ever, interestingly, in the presence of any disorder or fluctuations in the 
kinetic, structural or mechanical parameters, a window opens up be- 
tween the growth and collapse phases that allows for the coexistence of 
both states, which naturally characterizes the region of dynamic insta- 
bility in terms of experimentally measurable parameters. In addition 
to providing an explanation for the occurrence of dynamic instability 
in polymeric filaments, our study may be of possible significance for 
the evolvability of this trait. 

Mathematical model 

For simplicity of exposition, we consider a one-dimensional elastic fil- 
ament restricted to move in two dimensions and capable of attaching 
to or detaching from a rigid substrate with which the filament inter- 
acts via a series of springs (see Fig. [T^). The filament is assumed to 
be made of subunits of length a, and have a bending stiffness B and 
stretching stiffness E (where, as is usual B/Ea^ <C 1), while the linear 
springs that may connect it to the substrate are assumed to have a 
stiffness S and a maximum extension re- 

Given an initial state of the filament with say N segments that is 
attached to a substrate, we assume that its shape described by the posi- 
tions of the ends of its segments Vi{t) — {xi{t), yi{t)) evolves according 
to the equation 

df 
^r^--a^+f^W (1) 

corresponding to overdamped Langevin dynamics with rj the damping 
coefficient, and the elastic energy of the filament £ is given by 

1 



E^^d 



i=l 



2„„-.+i - r,| - a)Y -Bcos{e, - 0) + -S{r, - R0^ (2) 



where 9i are the angles between neighboring subunits, and Ri = {ai, 0) 
are the positions of the endpoints of the i—th spring on the substrate 
to which the filament is attached, with the proviso that S* = 0, if 
|ri — Ri| > Tc, i.e. the springs break when the extension equals or ex- 
ceeds the threshold r^. Here, the first term corresponds to the stretch- 
ing energy (which is negligibly small relative to the other terms), the 
second term is the bending energy associated with rotating one segment 
relative to its neighbor away from the natural rest angle between 
them, and the last term is the energy associated with the interactions 
with the substrate. The thermal fluctuations are characterized by fi(t), 
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with {ii(t)) — 0, and correlations 

(/7(i)//(i')> = 2T^kBT6,,6^sS{t - t'), (3) 

where f^ is the 7 component of the vector fi, T is the temperature, ks 
is the Bohzmann constant, and Sij is the usual Kronecker-delta oper- 
ator. The assumptions underlying the above dynamics are common in 
polymer physics: we ignore all long range hydrodynamic interactions, 
and furthermore assume that detailed balance and the fluctuation- 
dissipation theorem are valid at the scale of the filament. 

Subunits can attach to an end of the filament with a rate I/tq, 
where tq is a characteristic growth time. We note that for polar fil- 
aments, the rate of attachment on either end of the filament will be 
different; here we dispense with this difference for simplicity, but can 
easily account for the new effects that it will lead to, including tread- 
milling. Newly attached subunits start out collincar with the filament, 
but evolve to acquire an intrinsic curvature n = (p/ a (where (p is the 
intrinsic angle between segments) via a single step first-order kinetic 
process that models the conformation change associated with NTP 
hydrolysis. The associated evolution equation for the curvature reads 

-^ = — (ko-k), (4) 

at Tk 

where kq — (po/o. is the intrinsic curvature at equilibrium following 
hydrolysis and r^ is the characteristic timescale for the evolution of 
curvature. 

In addition to ambient thermal fluctuations, the structural, kinetic 
and mechanical properties of the filament in general could fluctuate in 
space due to various source of local heterogeneities in hydrolysis rates, 
or the adhesion to the substrate to name just two possibilities. This 
type of heterogeneity leads to quenched disorder (i.e. time indepen- 
dent disorder) and we will consider a particular manifestation of it as 
exemplified by a random distribution of the adhesion spring constants 
Si that leads to a random distribution of spring toughness along the 
filament. In particular we use the Gamma distribution given by 

p^^^ ^ — ji^nk) ' ^^^ 

where k is an integer and fi is the mean of the distribution (i.e. (S) = 
n). With this choice, changing k allows us to keep the mean of the dis- 
tribution constant and vary its standard deviation a = fJ,/k, allowing 
us to tune the degree of disorder and quantify it using the coefficient 
of variation Cy — a / ^ — 1/k. Recent experiments using varying con- 
centrations of non-hydrolyzable analogs of tubulin and their variants 
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in microtubules are consistent with this type of quenched disorder, as 
we discuss later. 

There are two natural dimensionless parameters in the problem: 
(i) a^ = Bk^ /Sri, a ratio of the linear energy density associated with 
straightening out a filament with a natural curvature k and the linear 
energy density associated with the adhering springs at their point of 
failure Sr^ and (ii) /3 — t^/tg a ratio of the time scales for hydroly- 
sis to subunit addition, in addition to the coefficient of variation that 
characterizes the disorder 1/fc. As we will see these two parameters 
characterize the phase space of possibilities for the stability and dy- 
namics of the filaments while the presence or absence of thermal and 
quenched disorder controls the appearance of a window of coexistence 
of the growing and shrinking phases where dynamic instability arises. 
In Table [U we list the main parameters employed in the model. 

To study the dynamics of Eq. [T]-Eq. 21 we prescribe an initial state 
of the filament with one end of the filament strongly attached to the 
substrate so that we need to track only the dynamics of the other end 
as it grows and shrinks. Our analysis starts with a consideration of the 
static problem before we progressively increase the level of complexity 
of the model. We start by considering the dynamics of the filament 
in the limit tq — ^ oo (i.e. no growth), in a fully hydrolyzed state at 
temperature T — and in absence of quenched disorder. We then sep- 
arately study the effects of quenched disorder and thermal fiuctuations. 
Finally, we consider the full model with growth and hydrolysis. In each 
case, we use a combination of analysis and numerical simulation; the 
latter is carried out using a conjugate gradient method (for the static 
case), and a fourth order adaptive step Runge-Kutta method in the 
limit r = and a simpler Euler scheme in the case T > (for the 
dynamic cases). 

Analysis 

Statics: unbinding transition 

We start our analysis of the mathematical model by focusing on the 
simplest of cases corresponding to the static solution of Eq. [T]at T = 
(i.e. without noise) which satisfies 

In this setting, we let tq be larger than all other time scales in the 
problem, so that the filament has a constant length. Furthermore 
we impose that all the subunits have a constant intrinsic curvature 
K — (jj/a. Then the system of equations Eq. [5] can be efficiently solved 
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using the method of conjugate gradients in a manner that is much 
faster than a dynamic simulation, an approach we use to simulate the 
model ioT B = 1,E = 10^, re = 1, a = 1, A'' = 128 where we note 
that B/Ea^ ^ 1 so that polymer stretching is unimportant. Then the 
behavior of the model only depends on the dimensionless ratio a, via 
the static subunit angle (j) = kg (which is constant here) and the lateral 
stiffness S. For small (j> the filament is completely attached; once (j> is 
increased past a critical value, the filament unbinds from the substrate 
starting from a free edge. The unbinding transition is abrupt as can 
be seen computing the steady-state value of the spatially averaged 
local bending angle ^ - Oj/N (or analogously the scaled curvature) as 
a function of the intrinsic bending angle 4>. At a critical value ac, itself 
a function of S, we see that the filaments unbinds globally, as shown 
in Fig. [TJj. 

In order to obtain an analytical understanding of the unbinding 
transition, we consider the continuum limit and restrict ourselves to 
small slopes and deformations, parametrizing the filament coordinates 
in terms of the horizontal component x: r^ — >■ {x,y{x)), where y{x) 
describes the filament shape. By neglecting stretching of the polymer, 
the quadratic energy functional can be written as 

]^{aB{y" - nf + Sy^)dx. (7) 

The Euler-Lagrange equations associated with minimizing £ lead to 
the equilibrium equation (Eq. ^ which now read 

By"" + 5y = 0, (8) 

The associated boundary conditions are 

y"{Q) = K By"'{0)^-Syi0)a/2, y{<x),y'{<x) ^ (9) 

Here the first boundary condition imposes the fact that the filament 
has an intrinsic natural curvature k, the second boundary condition 
accounts for the discrete nature of the adhering spings and arises by 
noting that the force from the lattice of springs can be rewritten as 
— Sa ^^ S{x — ia)y, while the last two simply state that a semi-infinite 
filament strongly attached at infinity is horizontal. Solving the bound- 
ary value problem Eq. (8-9) leads to a shape profile for the filament 
given by 

where q = [SjBYl^ I\f2 and characterizes a natural healing length Ih = 
1/(7 = \f2(B I S)^!^ . Thus over scales large compared to In from the free 
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edge, the filament remains straight when bound to the substrate. This 
linear solution agrees very well with the numerical solution as shown in 
Fig[TJ;. Not accounting for the discrete nature of the adhering springs 
is tantamount to setting a = 0, so that the solution Eq. [TU] reduces to 
y{x) = ^e~'^^ {cos{qx) — sin(qx)) which slightly underestimates the 
numerical solution, a result similar to that obtained in (|ll^ for a more 
complex cylindrical geometry. The analytical solution Eq. [TU] allows us 
to determine the condition for unbinding of the filament which would 
naturally start at the free edge when y{0) = r^ i.e. the spring at 
the edge has been stretched to its maximum length and would thus 
break. In the absence of fluctuations and disorder, this would cause 
the entire filament to unzip. If Kc is the critical natural curvature above 
which the filament will spontaneously unbind, we find that the critical 
condition reads y(0) = Kc/2g^(l + aS'o/(4J3g'^)) — r^ which can be 
recast in terms of the dimensionless parameter a = \J Bk^/ Sr"^ yielding 
a critical structural parameter a^ = (1 + p) which delineates the bound 
and unbound states. In the inset Fig. ^p, we show that the rescaled 
numerically obtained critical curvature Uc fits the simple theoretical 
expression above very well, and as the healing length lu/a — >■ oo, etc = 1 
as expected. To understand the dynamical process of unbinding, we 
now consider the effects of both deterministic and stochastic processes. 

Deterministic dynamics 

When a > «£, we numerically integrate the equations of motion Eq. 
(IlJ and follow the evolution of the length of the detached part of the 
filament, which is still assumed to be of uniform length over the time 
scale of observation. The results, shown in Fig. [2^, suggest that the 
detached length grows diffusively in time. This behavior can be ex- 
plained by a simple scaling argument that balances the viscous dissi- 
pation rate Py ~ r]{l/t)'^l, where I is the length of the detached filament 
with the elastic power that drives unbinding Pe ~ {Bk^ — Sr'^)l/t = 
B{k^ - {KcY)l/t, and yields / ~ (^(k^ - {K^f)t/f]Y/'^. The slowing 
down of unbinding with time arises because the ever lengthening un- 
bound part takes a longer and longer time to move through the viscous 
environment; eventually once the unbound part has formed a circular 
ring, this diffusive behavior will likely be replaced by linear Stokesian 
dynamics, although we do not reach this limit in our simulations. In 
contrast, if the sub units break off as soon as they unbind from the 
substrate, viscous dissipation is localized to a region near the dynamic 
detachment zone, and the viscous dissipation rate Py ~ rj{l/t)'^a so 
that now I ^ B{k^ — {Kc)'^)t/ar]. To verify this relation, we use numer- 
ical simulations where we remove the subunit when it detaches from 
the substrate; the results shown in Fig. [21d confirm that we indeed 
capture this Stokesian limit as well. 
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Dynamics with quenched disorder 

When we introduce quenched disorder into the structural parameter ac 
via the dependence of the substrate stiffness S through its coefficient 
of variation C„, unbinding occurs stochasticahy. In Fig. [3^, we show 
the resuhs of simulations for the steady-state value of the average lo- 
cal bending angle {J2n ^j/-^) a-s a function of the intrinsic bond angle 
(j) = an, for different values of the coefBcient of variation Cy, keeping 
the average stiffness constant (in this example /j, = 0.05). Since every 
calculation starts with a particular realization of the quenched disor- 
der, our results are shown as averages over different realizations of the 
disorder. 

We note that disorder causes the unbinding of the filament to occur 
for values of (j) that are larger than when disorder is absent {Cy = 0). 
The seemingly counter-intuitive result that disorder makes the system 
stronger is immediately rationalized once we realize that this is only 
true because unbinding is always ruled by the strongest region, unlike 
material fracture or failure that is controlled by the weakest bond (J20h . 
For N independent random variables Si, distributed according to Eq.[5l 
the probability that the maximum is less than S is equivalent to the 
probability that all the values 5*^ are less than S. Hence, the cumulative 
distribution of the maximum over N values is given by 

PNiS)^PiSf (11) 

where P{S) = J^ p{x) dx is the cumulative distribution associated to 
p{x). We can estimate the unbinding threshold considering an equiv- 
alent filament attached with springs whose stiffness is given by the 
average Smax of the distribution in Eq. 1111 The result, in dimension- 
less terms, reads as 



Qfc ~ (1 + a/lh)y/Smax/fJ- (12) 

Our dynamical simulations allow us to deduce this relationship numer- 
ically, and as shown in the inset of Fig. [3ji, there is good agreement 
with the simple theory. It is also worth pointing out, that in the ther- 
modynamic limit of large N, the asymptotic limit of Eq. [TT] is given 
by the Gumbel distribution: 

Pn{S) ~ exp^A^e-'^^/^] (13) 

Then the average value of the stiffness is given by {S,nax) — /^(t + 
log N)/k, where 7 ~ 0.57 is the Euler constant, or in dimensionless 
terms, a^ - (1 + a/lh)\/il + logN)/k. 

Dynamics with thermal fluctuations 

In the presence of thermal fiuctuations but without any quenched dis- 
order, the unbinding of the filament can be activated by noise even 
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though the filament is nominally stable from a deterministic perspec- 
tive. This is because, at sufficiently low temperatures, one can reason- 
ably assume that unbinding always starts from the edge of the filament 
and proceeds by breaking the bonds sequentially Q. Then suggests that 



the typical unbinding velocity follows the Arrhenius law (J2lh 



where the energy barrier is given by AE ~ ^('^''c ^ Bn^jlh ^ ac — a. 
Our numerical simulations of Eq. [T]in a system with a = O.bac and 
different values of the temperature T confirms this picture, as shown 
in Fig. [Id. 

Kinetics with thermal and quenched disorder 

We finally tie all the elements of the minimal model together to study 
the kinetics of an adsorbed polymorphic filament that can grow and 
change conformation as characterized by the two natural control pa- 
rameters in the problem: /3 = t^/tg the ratio of the rate of confor- 
mation change (NTP hydrolysis) to the rate of addition of subunits 
and the maximum value reached by the structural/mechanical con- 
trol parameter ao (where the subscript now denotes the equilibrium 
value of the parameter, which is also its maximum, in contrast with 
its dynamical value a). A natural link between the two dimension- 
less parameters in the problem can be seen by rewriting Eq. 2] as a as 
da/dt — Ti^{ao — a), where ao = KovB/rc-^/JI and n is the average of 
the spring stiffness S. 

In Fig.|4]we show the evolution of the scaled filament length L/a for 
different value of ao and /3, with springs of random stiffness, distributed 
according to Eq.[S]with k — -i and /i — 0.05. In our simulations, we see 
that for constant f3 the filament grows for small ao and collapses for 
large ao (Fig. Hji). Similarly, for constant ao, the filament grows for 
small (3 and collapses for large /3 (Fig. [JJd). However, for intermediate 
ranges of ao,/3, we observe a coexistence of the growth and collapse 
regimes, punctuated by intermittent rescue events. In this case, we see 
that the rescue events often occur at the same location, consistent with 
our earlier analysis which shows that a single tough bond is sufficient to 
stop the fracture from propagating. However, once a new tough bond 
is formed further along the growing filament, it serves to arrest the 
next catastrophe. The heterogeneous nature of the observed dynamic 
instability in our minimal model is related to the presence of quenched 
disorder in the mechanical parameter ao, since the kinetics are assumed 



^internal breakages are of course possible, but they result in self- limiting catastrophes 
since any internal bubble of detachment will lead to a weakening of driving stress. 
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to be completely deterministic. These results are consistent with recent 
experiments that show that rescue events are directly correlated with 
the presence of GTP-bound subunits in a microtubule (|l9i ) , which serve 
to arrest catastrophes and rescue them repeatedly at the same location. 

If instead, we use a structurally homogeneous system with no quenched 
disorder but with stochastic hydrolysis events, then Eq. |4] describes 
the process only on average. Allowing each subunit to be hydrolyzed 
stochastically with rate rn = ^/tk,, switching its intrinsic curvature 
from K = to K = Ko (or equivalently from a = to a = ao) leads 
to a crossover from growth to collapse as (3 is increase, as shown in 
Fig-Stj but with one important difference - in the absence of quenched 
structural disorder, catastrophes are much more likely to go all the 
way to collapse. Of course varying the the type of stochasticity in the 
hydrolysis will likely lead to different types of intermittency. 

In Fig. \5\ we show a phase diagram showing growth, collapse and 
dynamic instability as a function of ao and (3 for the case of structural 
disorder and deterministic hydrolysis. The boundary between growth 
and collapse in the fully deterministic limit can be understood using 
simple considerations by noting that the filament will collapse when 
its curvature reaches the unbinding point before a sufficient amount 
of stabilizing subunits are attached to its end. The unbinding point is 
reached in a time tc given by a{tc) = ac and we can thus write the 
condition for collapse as tc = utq, where n is a constant. Inserting 
this condition in the solution of Eq. 21 we obtain an estimate for the 
phase boundary as 

that is in qualitative agreement with the numerical results. Adding 
quenched disorder opens up a lenitcular region along this boundary 
where co-existence of the growth and collapsed phases leads to dy- 
namic instability. In our analysis, we have assumed that the filament 
dynamics is much faster than the growth kinetics and the evolution 
of the structural parameter a. If we relax this constraint, our model 



could also show a regime, defined "third state" in Ref. (|ll[ ) where 
catastrophes are interrupted by the attachment of new segments; this 
will be the subject of future work. 

Discussion 

Inspired by the occurrence of metastability of a growing and shrinking 
phase in a variety of polymeric assemblies, we have studied a minimal 
model of a kinetically polymorphic elastic filament adsorbed on a soft 
substrate that exhibits a window of dynamic instability. In particular, 
we find that the existence of this window depends on the existence of 
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disorder in a phase space that controls the growth and shrinkage of 
the filament. We characterize the range of behavior of the filament in 
terms of two parameters, a structural-mechanical parameter that char- 
acterizes the stored energy in the residually strained assembly, and a 
kinetic parameter that characterizes the relative rate of growth to the 
rate of internal conformation change (associated with NTP hydrolysis). 
In the absence of any disorder, we find that the criterion for filament 
to remain bound or not is characterized by a structural parameter that 
measures the relative strength of intra-filament to filament-substrate 
bonds. This parameter has a natural extension to multi-filament poly- 
meric assemblies such as microtubules, actin and their homologs in 
terms of intra-filament and inter-filament interactions which are known 
to play an important role in the mechanics and dynamics of these ob- 
jects. Furthermore, we find that if the filament unbinds, it does so 
via either diffusive or Stokesian dynamics, depending on whether the 
filament does not or does disassemble as it unbinds. 

The presence of disorder, either quenched or thermal, changes the 
unbinding picture qualitatively. Quenched disorder in the structural- 
mechanical parameter effectively increases the unbinding strength, a 
result that initially seems counter-intuitive given that disorder nor- 
mally makes a sample weaker and thus more susceptible to fracture or 
peeling. Here, in this one-dimensional system, the effect is opposite 
due to the fact that stability is ruled by the strongest link. Thermal 
fluctuations on the other hand lead to subcritical unbinding, with a 
creep velocity that follows a simple Arrhenius law as expected from 
studies of other similar systems (|2ll ) . Although our treatment of the 
kinetic parameter is rather simple, it is able to capture the competi- 
tion between the rate of subunit addition and internal conformation 
change, a process implicated in the existence of a GTP cap; here we 
show that this leads to natural undulatory shape near the edge similar 



to that proposed earlier for more complex models (jlll ) . When we con- 
sider the kinetics of conformation change in addition, we find that a 
window of dynamic instability opens up at the boundary that separates 
unlimited growth from collapse in the presence of quenched disorder. 
In this region, the growth and shrinkage phases co-exist and our phase 
diagram for its existence points directly to the quantitative conditions 
for the existence of dynamic instability for our simple model. This pair 
of dimensionless measures which characterize structural stability and 
the kinetics of assembly and conformation change in the presence of 
randomness may be generalized to more complex geometries such as 
sheets, tubes, helices and beyond and raise an interesting question with 
evolutionary implications: what combination of geometry, mechanics 
and chemistry can lead to the conditions for dynamic instability, and 
thus the flexibility to build structures that are functional, and yet not 
permanent ? Understanding the chemical constraints on the ranges of 
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ao and /3 is clearly an important next step, and an artificial approach to 
this might well involve using a combination of known non-hydrolyzable 
analogs of GTP (NTP) along with the ability to track the locations 



(|22f ) associated with quenched structural disorder. 
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Figure 1: Filament conformation - statics, a) Schematic representation of 
the system consisting of an elastic filament with a non-zero natural curvature 
adhered to a substrate via a series of springs. Here the filament is akin 
to a protofilament while the substrate represents the "bath" of the other 
filaments, b) Static phase diagram showing the unbinding threshold in terms 
of the average local bending angle as a function of the intrinsic angle (f) for 
different values of the stiffness S. The inset shows the critical value of the 
dimensionless parameter a, that characterizes the ratio of the filament and 
substrate energy, as a function of the ratio between the healing length Ih and 
the discretization step a. The solid line is the theoretical result, in excellent 
agreement with the simulations, while the dashed line is the asymptotic 
result for a — )■ 0. c) The characteristic shape of the attached filament in 
the neighborhood of the edge shows the healing length characterizing the 
balance between filament and substrate deformation. The simulation results 
are compared with our linearized continuum theory. 
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Figure 2: Filament dynamics for different values of the intrinsic curvature 
K and spring stiffness S, corresponding to the unbound phase a > Oc- The 
filament unbinds dynamically - there are two possible mechanisms that de- 
termine the kinetics of the process, (a) Diffusive unbinding occurs when 
the filament detaches partially but remains intact. In this case, the figure 
shows that the detached length increases as l/a = V*(t/r])'^''^. The inset 
shows that the coefficient V* , for different values of S, follows a single linear 
function when plotted as a function of (k^ — KcY''^- (b) Stokesian unbinding 
occurs when the subunits break off from the filament once they are detached. 
The Figure shows that the detachmed length grows as l/a = V{t/rj)/a. The 
inset shows that the coefficient V scales as is a linear function of (k^ — k^) 
as expected for small to moderate values of the parameters. 



Figure 3: Filament dynamics in the presence of quenched disorder and ther- 
mal fiuctuations. (a) When quenched disorder is introduced in the spring 
stiffness, the unbinding threshold in Figure 1 is broadened due to the statis- 
tics of extremes as characterized in terms of the coefficient of variation 
Cv = o"/)U, the ratio of the standard deviation to the mean of the spring 
stiffness S. Inset shows the critical value of the dimensionless parameter 
a = \/Bk/ ^Jjivc as a function of the coefficient of variation of the disor- 
der. The error bar represents the variance of oic rather than the error on 
the mean. Both the mean and the variance of Oc are seen to increase with 
increasing C„ as predicted by extreme value statistics (line), (b) Due to 
thermal fluctuations, the filament can unbind also for a < oic- Here, we 
show the decrease of the filament length for a = 0.5ac. The inset shows 
that the velocity in the initial stage satisfies the Arrhenius law. 



Figure 4: Filament growth and collapse kinetics, a. Filament length as 
a function of time for a range of ao the mechanical control parameter, b. 
Filament length as a function of time for a range of /3 = t^Jtg^ the ratio 
of the time constant for the intrinsic curvature to equilibrate relative to the 
time constant for subunit addition, the biochemical control parameter, c. 
Filament length when the hydrolysis is random, also shows the transition 
between growth and catastrophe, but with one important difference from 
the case when the toughness parameter a is random; shrinkage goes all the 
way to collapse in the absence of any structural mechanical inhomogeneities. 
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Figure 5: A phase diagram in a— /3 space summarizing the regimes of growth, 
collapse and dynamic instability, a window that appears in the presence of 
disorder either in the structural toughness of the filament-substrate bond, or 
the hydrolysis. The dashed line is the theoretical prediction corresponding 
to Eq. dSl 



Parameter 


Description 


Elastic constants: 


B 

E 
S 


filament bending stiffness 
filament stretching stiffness 
spring stretching stiffness 


Characteristic lengthscales: 



re 

a 

L 



h = 1/g = V2(B/Sf'^ 



spring maximum extension 

segment length 

length of attached portion of the filament 

length of detached portion of the filament 

healing length 



Angles and curvatures: 



aK 



rip 



intrinsic curvature 

intrinsic angle 

local filament angle 

maximum value of k when k is evolving 

unbinding threshold 



Dynamics: 






filament damping coefficient 
timescale of curvature evolution 
timescale of filament growth 



Disorder: 



a 

C^ = a 111 



average spring stretching stiffness S 
standard deviation of stiffness distribution 
coefficient of variation 
maximum value of S along the filament 



Control parameters: 



a = yBn/\/Src 
P = t^/tg 



mechanical control parameter 
maximum value of a when k is evolving 
kinetic control parameter 



Table 1: A description of the parameters of the model. 
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